Particle propagation in cosmological backgrounds 



Daniel Arteaga 

Departament de Fi'sica Fonamental. Facultat de Ffsica. Av. Diagonal 647, 08028 
Barcelona (Spain) 

E-mail: darteagaOub . edu 

(N 

O I Abstract. We study the quantum propagation of particles in cosmological 

' backgrounds, by considering a doublet of massive scalar fields propagating in an 

expanding universe, possibly filled with radiation. We focus on the dissipative effects 
related to the expansion rate. At first order, we recover the expected result that the 
decay rate is determined by the local temperature. Beyond linear order, the decay 
O -1 rate has an additional contribution governed by the expansion parameter. This latter 

| ' contribution is present even for stable particles in the vacuum. Finally, we analyze the 

long time behaviour of the propagator and briefly discuss applications to the trans- 
Planckian question. 
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In this contribution we study the quantum propagation of particles in a cosmological 
\ background. We are particularly interested in understanding the dissipative phenomena 
related to the time dependence of the metric. To this end, we analyze the propagator 
of a massive particle which interacts with a massless radiation field in an expanding 
universe. Several points must be considered. 

First, we are dealing with an interacting field theory in a curved spacetime. In 
this situation, the asymptotic in and out vacua generally do not coincide. Being 
interested in expectation values, rather than in in-out matrix elements, we adopt the 
Keldysh-Schwinger formalism, or Closed Time Path (CTP) method [1-3] in curved 
spacetime [4-7]. 

Second, as it is well known, in a curved spacetime there is no single definition for 
the vacuum nor for the concept of particle. We face this issue by working within the 
adiabatic approximation [8] : the massive particles will have their Compton wavelengths 
much smaller than the typical curvature radius of the universe (in our case, the Hubble 
radius) . 

Third, as explained in [9], in theories such as QED or perturbative quantum gravity, 
dissipative effects appear only at two loops, because the one-loop diagrams which could 
lead to dissipation vanish on the mass shell. Here, in order to keep the calculations 
simple, we have chosen a simple, yet physically meaningful, model which exhibits 
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dissipation at one loop. We expect the behaviour of QED or perturbative quantum 
gravity to be similar at two loops. 

In this contribution we compute the retarded self-energy of the lightest field in a 
massive doublet which propagates in a thermal bath of massless particles in an expanding 
universe, and from it we extract the decay rate. Notice that in the Minkowski vacuum 
the excitations of the lighter field are stable, hence the decay rate is zero. We concentrate 
on the physical insights and summarize the main results. A more detailed account will 
be given in separate publications [10, 11] 

The contribution is organized as follows. In section 2 we introduce the model and 
motivate the use of the adiabatic approximation for the massive fields. In section 3 
we present the results for the imaginary part of the self-energy and the decay rate. In 
section 4 we study the time evolution of the interacting propagator. Finally, in section 
5 we summarize the main points of the contribution and discuss its relevance to the 
trans-Planckian question. We use a system of natural units with h — c — 1, and the 
metric has the signature (— , +, +, +). 



2. The model 



We consider spatially isotropic and homogeneous Friedmann-Lemaitre-Robertson- 
Walker models with flat spatial sections: 

ds 2 = -dt 2 + a 2 (t)dx 2 . (1) 

The particle model is the following: two massive fields (f> m , and (pM, interacting with a 
massless field, x> via a trilinear coupling. The total action is S = S m + Sm + S x + S- mt , 
where each term is given by 

S m = \f dt d 3 xa 3 (t) ((<9 t m ) 2 - -L^(d x( f> m ) 2 - m 2 0^ , (2a) 

S M = \J dtd 3 xa 3 (t) (W M ) 2 - -Ljid^M) 2 - M 2 ^ , (26) 

S x = \j dtd 3 xa 3 (t) (W) 2 - -^y(^xX) 2 -m)x 2 ) , (2c) 

S int = gM dtd 3 xa 3 (t)0 m 0MX, (2c?) 



with R(t) being the Ricci scalar. We assume that the massless field is conformally 
coupled to gravity, so that £ = 1/6. It is useful to work with rescaled massive fields 
defined by 0(t,x) := [-g(t, *)} 1/4 <f)(t, x) = a 3 / 2 (t)0(t, x). 

We consider the two massive fields having large masses but with a small mass 
difference Am := M — m M. As shown in [12], the model can be interpreted as 
a field-theory description of a relativistic two-level atom (of mass m and energy gap 
Am) interacting with a scalar radiation field x- The radiation field x is assumed to be 
at some conformal temperature 9 (which can eventually be zero). The corresponding 
physical temperature, as well as the Hubble rate H(t) := a(t)/a(t), are chosen to be 
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much smaller than the masses of the fields. These restrictions ensure that the number of 
massive particles is strictly conserved. The non-trivial dynamics concerns the transitions 
between the two massive fields accompanied by emission and absorption of massless 
quanta. 

In a curved spacetimes it is not a trivial task to compute even the free field vacuum 
propagators. For massless conformally coupled fields there is a natural vacuum state, 
the conformal vacuum. Propagators in this vacuum, when expressed in conformal time, 
essentially correspond to the flat spacetime propagators. For the massive fields, rather 
than attempting to find the exact free propagator, we will exploit the fact that their 
Compton wavelengths is much smaller than the Hubble length H~ x . In this regime, the 
adiabatic (WKB) approximation is valid and explicit expressions for the free propagators 
can be computed [8,10] — see for instance ( fT9l ). 

3. The self-energy and decay rates 

In this section we consider the interacting retarded Green function 

g r (m';p) :=e(t-t'){^ mp (t)^ mp (t')}) (3) 

within the adiabatic approximation. It is related to the retarded self-energy £ R via [12] 
G R (M';p)=<4 0) (M';p) 

-i j ds ds' y/-g(s) V-9(s')G^ (t, s; p)E R (s, s'; p)G R (s', t'; p) (4) 

where G^(t,t';p) is the free retarded propagator. In terms of the rescaled fields, 
4>(t; p) = a 3//2 0(t; p), the above relation becomes 

G R (t,t';p)=G { °\t,t';p)-t J dads' S ;p)E R ( S , s';p)G K (s',t';p). (5) 

We assume that the massive fields are in the adiabatic vacuum, and that the 
massless field \ is in a thermal state, characterized by a fixed conformal temperature 
9. We will compute the imaginary part of the one-loop self energy to order g 2 in 
the adiabatic approximation, evaluated at the mass shell. It will be evaluated in a a 
frequency representation around the average time coordinate T = (ti+t 2 )/2, by Fourier- 
transforming with respect to the difference coordinate A = — t 2 , which amounts to a 
local frequency representation (it is further analyzed in next section). As for the spatial 
part, we work in the momentum representation to exploit conservation of the conformal 
momentum. 

3.1. Linear approximation to the scale factor 

As a first step, we approximate the evolution of the scale factor by a linear expansion: 

a(t) na(T)[l + H(T)(t-T)}. (6) 

This approximation for the scale factor is appropriate when considering physical 
temperatures which are much larger than the expansion rate (but still much smaller than 
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the fields masses). On the mass shell, thermal corrections will dominate over curvature 
corrections, since the thermal energy scale is much larger than the curvature energy 
scale. Therefore, we expect the on-shell self-energy to be governed by the thermal bath 
at the instantaneous physical temperature at each moment of the expansion, 6/a(T). 

The explicit calculation [10] confirms that the imaginary part of the on-shell self- 
energy is given by that of a thermal bath in Minkowski at a physical temperature 9/a(T). 
In the limit in which the atoms are at rest this result is [10, 12] 

Q 2 

Im£ R (m,T;0) = -^-MAmn e/a{T) (Am), (7) 
where n e / a ( T )(Am) is the Bose-Einstein function: 

n e/a(T) (Am) := eAma ^ T)/e _^ (8) 

As in Minkowski spacetime, the self-energy corresponds to a decay rate, 

1 Q 2 
r = Im E R (m, T; 0) = §- Am n 0/a(T) (Am) , (9) 

which amounts to the probability per unit time for the lightest state to absorb a massless 
particle from the thermal bath. 



3.2. Beyond linear order: vacuum effects 

When the expansion rate of the universe is of the order of the temperature or larger, 
vacuum effects become relevant. Energy conservation does not hold for energy scales of 
the order of the expansion rate, and therefore we expect new channels for the particle 
decay which will contribute to the imaginary part of the self-energy. 

In order to study the vacuum effects we need to choose a explicit model for the 
evolution of the scale factor. For instance, in the case of de Sitter, 

a(t) = a(T) e H ^ T \ (10) 

the vacuum contribution to the imaginary part of the retarded self-energy given by [11] 

Q 2 

Im£ R (m,T;0) = -f-MAron H/(2ff) (Am), (11) 

o7T 

which coicides with the self-energy in a Minkowski thermal bath at a temperature 
H/(2ir). The result is not unexpected since the effective de Sitter temperature [13] 
is recovered. The corresponding decay rate 

1 - a 2 
T = imE R (m,T;0) = f-Amn ff/(2ff) Am , (12) 

771 07T 

amounts for the probability per unit time for the lightest field to emit a massless particle. 
Energy conservation forbids this process in Minkowski spacetime, but this restriction 
does not apply in an expanding universe. 
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4. Retarded propagator and self-energy in cosmology 

In expanding universes the propagators are no longer time-translation invariant. We 
can nevertheless always express the propagator in a frequency representation, 

G R (cu, T; p) := J dA e^ A G R (T + A/2, T - A/2; p) . (13) 

For short time differences as compared to the inverse expansion rate, i.e., \t — 1'\ <C H^ 1 , 
([5]) can be diagonalized: 

_ — 2 

C R (a;,r;p)= (14) 
[-«GW(w,T;p)] 1 + S R (^,T;p) 

Fourier-transforming again we get the short-time behavior: 

<5r(m';p) = 1 fL-,sm[R P (T)(t-t)) e-^™-^e(t-t'). (15) 

with 

R 2 p (T) := E 2 p (T) + Re E R (E p , T; p) := m 2 + + Re S R (£ p , T; p) (16) 

and 

r p (T) := -^-L^ImS R (E p ,T;p). (17) 

Therefore one recovers the usual interpretation, in which the real part of the self-energy 
corresponds to the energy shift, and in which the imaginary part corresponds to the 
decay rate. Notice that both quantities depend in general on time. 

One may also be interested in considering large time lapses, and in this case the 
frequency representation of the propagator around the average time does not make sense. 
Lifting the short-time requirement, and only imposing the adiabatic approximation, the 
following expression for the evolution of the retarded propagator is found [10]: 

G R ( tl ,t 2 ;p) = ~\ sin ( r df R p (t')) e-fiWWeitt-tt). (18) 

V -ftp 1*1 J -ftp (*2 J \Jt 2 J 

Notice that the long-time evolution of the propagator can be expressed in terms of 
integrals of quantities evaluated in the local frequency representation. Two time scales 
are clearely separated: the interaction timescale, in which the interaction process take 
place and in which the self-energy is evaluated, and the evolution timescale, which 
can be much longer and during which the propagators deviate significantly from the 
corresponding Minkowski expression. Equation (fT8|) can be derived in a very similar 
way as the well-known adiabatic approximation for the free retarded propagator [8]: 

<3g°(ti,t2;p) = /F = n sin ( f 1 dt' E p (t')) 9{t x -t 2 ). (19) 
V-^p^iJ-^pv^J \Jt 2 J 
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5. Summary and discussion 

The goal of this contribution is to analyze the quantum effects in the propagation of 
interacting fields in a cosmological background. This issue may play an important role in 
justifying the non-trivial dispersion relations which have been used when addressing the 
trans-Planckian question in the context of black holes [14-17] and cosmology [18-21]. 
Interactions could indeed significantly modify the field propagation when approaching 
the event horizon of a black hole [22-25] or at primordial stages of inflation [9]. 

In our model, the masses of the fields were assumed to be much larger than the 
expansion rate of the universe. This was a key assumption, because it allowed to 
introduce the adiabatic (WKB) approximation, which not only makes the problem 
solvable, but also allows having a well-defined particle concept even in absence of 
asymptotic regimes. Within this approximation, the time-evolution of the interacting 
propagators can be computed from the integral of the retarded self-energy, evaluated 
on-shell in a frequency representation around the mid time. 

The imaginary part of the self-energy determines the decay of the retarded 
propagator, and hence it is an expression of the dissipative properties. For temperatures 
higher than the expansion parameter the decay of the propagator is determined by the 
local temperature at each moment of expansion. For lower temperatures, the decay of 
the propagator is driven by the expansion rate of the universe. This second contribution, 
which is present even in the vacuum, can be interpreted as being a consequence of the 
absence of energy conservation at those energy scales comparable to the expansion rate. 

The decay rate, derived from the imaginary part of the self-energy, has a secular 
character. Even small decay rates could thus give an important effect when integrated 
over large periods of time. The exact significance of the generically dissipative properties 
of the propagator will be further analyzed elsewhere [11]. 
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